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Abstract 

We study the Hankel transforms of sequences related to the central coefficients of a 
family of Pascal-like triangles. The mechanism of Riordan arrays is used to elucidate 
the structure of these transforms. 
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1 Introduction 

This note concerns the characterization of the Hankel transfoms of the central coefficients 
T(2n, n, r) of a family of Pascal-like triangles that are parameterised by an integer r. Specif- 
ically, we define a family of number triangles with general term T(n, k, r) by 



*M.r,-gg)( 



k\ (n — k 



For instance, r = 1 gives Pascal's triangle A007318 , while r = 2 gives the triangle of Delannoy 
numbers [1], A008288 . 

Proposition 1. The Hankel transform of the sequence a(n,r) = T(2n,n,r) is given by 

fn + l\ 

2 n A 2 ). 

Proof. We proceed as in [9] and [7] by means of the LDL T decomposition of the Hankel 
matrix H(r) of T(2n,n,r). We take the example of r = 2. In this case, 
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Then 
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Hence the Hankel transform of T(2n, n, 2) is equal to the sequence with general term 

n 

Y[(2.2 k - k ) = 2 n 2( r 

fc=0 

L{2) is in fact the Riordan array 



1 — 3a; — \/l — 6 a; + 



x- 



or 



yl — 6x + x 2 
, 1 - 2x 2 



Ax 



x 



v 1 + 3a; + 2a; 2 ' 1 + 3x + 2a; 2 ' 
In general, we can show that H(r) = L(r)D(r)L(r) T where L(r) is the Riordan array 

1 1 - (r + l)x - y/l - 2(r + l)x + (r - l) 2 x 2 



(- 



2rx 



v /l-2(r + l)a; + (r-l) 2 x 2 ' 

and D(r) is the diagonal matrix with n-th term 2.r n — n . Hence the Hankel transform of 
T(2n, n, r) is given by 



f[{2.r k - k ) = 2 n r( n 2 1 ). 



fc=0 



□ 



We note that the Riordan array L(r) 

1 1 - (r + l)x - y/1 - 2(r + l)x + (r - l) 2 a; 2 , 



V 1 - 2(r + l)x + (r - l) 2 x 2 ' 
is the inverse of the Riordan array 

1 — rx 2 



2rx 



x 



1 + (r + l)x + rx 2 ' 1 + (r + l)x + rx 2 ' 
Its general term is given by 
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n\ ( n 
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(r + 1) 
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Its fc-th column has exponential generating function given by 

e( r+1 ^'4(2 v ^x)/v^ fe . 
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Corollary 2. The sequences with e.g.f. 7o(2yrx) have Hankel transforms given by 2 n r^ 2 /. 
Proof. By [1] or otherwise, we know that the sequences T(2n, n, r) have e.g.f. 

e {r+1)x I {2^x). 

By the above proposition and the binomial invariance property of the Hankel transform [6], 
B _r-1 T(2n, n, r) has the desired Hankel transform. But B~ r ~ 1 T(2n, n, r) has e.g.f. given by 



e- (r+1)x e (r+1)x /o(2 v ^x) = I (2y/r~x) 



□ 



2 Hankel transform of generalized Catalan numbers 

Following [1], we denote by c(n;r) the sequence of numbers 

c(n; r) = T(2n, n, r) — T(2n, n + 1, r). 
For instance, c(n; 1) = c(n), the sequence of Catalan numbers A000108 . We have 
Proposition 3. T/ie Hankel transform ofc(n;r) is rv 2 /. 

Proof. Again, we use the LDL T decomposition of the associated Hankel matrices. For 
instance, when r = 3, we obtain 
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Then 
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Hence the Hankel transform of c(n; 3) is 

n 



3^ 2 ) 



k=0 
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In this case, L(3) is the Riordan array 



V l + 3:?;' 1 +4x + 3x 2J 
In general, we can show that H(r) = L(r)D(r)L(r) T where 

ir r \ _ f 1 ^ \-i 

1 + rx ' 1 + (r + + rx 2 

and .D(r) has n-th term r n . Hence the Hankel transform of c(n; r) is given by 

n 

J r fc _ r ( 2 ) . 

k=0 

□ 

We finish this section with some notes concerning production matrices as found, for 
instance, in [4]. It is well known that the production matrix P(l) for the Catalan numbers 
C{n) = c(n, 1) is given by 

/ 1 ... \ 
1 1 ... 
1 1 1 ... 



pa: 



V 



/ 



Following [4], we can associate a Riordan array Ap(1) to -P(l) as follows. The second column 
of P has generating function Solving the equation 



u 



xu 



we obtain u(x) = — ^- — - = c(x). Since the first column is all O's, this means that Ap(1) is 
the Riordan array (l,xc(x)). This is the inverse of — x)). We have 



Ap(l) 
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10 

110 

2 2 1 
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Multiplying on the right by B, the binomial matrix, we obtain 

\ 



Ap(l)B 



( 1 
110 
2 3 10 
5 9 5 1 

V ; ; • ; 



L{1) 
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which is the Riordan array 



xc(x) 



1 — x' v l + x'l + 2x + x 2 

Similarly the production matrix for the c(n; 2), or the large Schroeder numbers, is given 



by 



P{2) 



/ 2 ... \ 
1 2 ... 
1 1 2 ... 
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• 7 



Here, the generating function for the second column is \zr- Now solving 

2 — xu 



u 



1 — xu 



which gives u = 1+x ^ 6rz -±^. Hence in this case, A P (2) is the Riordan array (1, — ^\ 6q ±g^' 
That is, 



A P {2) 
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2 4 

6 8 8 



V 



The row sums of this matrix are 1,2,6,22,90, 
right by the binomial matrix B, we obtain 



(1 fS_E))-i. 

2 — x 
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as expected. Multiplying Ap{2) on the 



A P (2)B 
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which is the array 



1—x — y/1 — Qx + x 2 1 — 3a; — y/1 — 6x + 



x- 



2x 



2x 
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Finally 
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which is 



or 



1 — x — yl — 6x + x 2 1 — 3x — yl — 6x + x 2 



2x 



Ax 



L(2) 



1 



x 



1 + 2x' 1 + 3x + 2x 2/ 
We can generalize these results to give the following proposition. 

Proposition 4. The production matrix for the generalized Catalan sequence c(n; r) is given 
by 

( r ... \ 
1 r ... 
1 1 r ... 



P(r) 



■ J 



The associated matrix Ap(r) is given by 



x(l-x) w 1 + (r - l)x - y/l - 2(r + l)x + (r - l) 2 x 2 

r — (r — lja; 2 

T/ie matrix L{r) in the decomposition L{r)D{r)L{r) T of the Hankel matrix H(r) for c(n; r) 
which is equal to Ap(r)B(l,x/r), is given by 



L(r) 
We have 



1 - (r - l)x - v 7 ! - 2(r + l)x + (r - l) 2 x 2 1 - (r + l)x - ^1 - 2(r + l)x + (r - l) 2 



x' 



2x 



2rx 



)• 
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1 + rx 1 + (r + l)x + rx 2 ' 



We note that the elements of L{r) 1 are in fact the coefficients of the orthogonal poly- 
nomials associated to H(r). 

Proposition 5. The elements of the rows of the Riordan array (j^, 1+ ( r+ f) x+ra .2 ) ar ^ the 
coefficients of the orthogonal polynomials associated to the Hankel matrix determined by the 
generalized Catalan numbers c(n;r). 



3 Hankel transform of the sum of consecutive general- 
ized Catalan numbers 

We now look at the Hankel transform of the sum of two consecutive generalized Catalan 
numbers. That is, we study the Hankel transform of c(n; r) + c(n + 1; r). For the case r = 1 
(the ordinary Catalan numbers) this was dealt with in [3] , while the general case was studied 
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in [8]. We use the methods developed above to gain greater insight. We start with the case 
r = 1. For this, the Hankel matrix for c(n) + c{n + 1) is given by 

/ 2 3 7 19 ... \ 

3 7 19 56 

H= 7 19 56 174 

19 56 174 561 

V ; ; 

Proceeding to the LDL T decomposition, we get 

H = LDL T 
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This indicates that the Hankel transform of c(n) 
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F(2A; + 3) 



/ V : : 

+ c(n + 1) is given by 
F(2n + 3). 



F(2k + 1) 

This is in agreement with [3]. We note that in this case, L -1 takes the form 
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where we see the sequences F(2n + 1) and (— l) n F(2n + 2) in evidence. 
Now looking at the case r = 2, we get 

/ 3 8 
8 28 
28 112 
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Proceeding to the LDL T decomposition, we obtain 
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Thus the Hankel transform of c(n; 2) + c(n + 1; 2) is 3, 20, 272, 7424 . . .. This is in agreement 
with [8]. We note that different factorizations of L~ l can lead to different formulas for h n (2), 
the Hankel transform of c(n; 2) + c(n + 1; 2). For instance, we can show that 

/ 1 ... \ / 1 ... \ / 1 
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We note that the diagonal elements of the last matrix correspond to the sequence a(n) of 
terms 1, 3, 5, 17, 29, 99, . . . with generating function 

1 + 3x — x 2 — x 3 
1 — 6a; 2 + x 4 

This is A079496 . It is the interleaving of bisections of the Pell numbers A000129 and their 
associated numbers A001333. We have 



a[n) 
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fc=0 



n + 1 
2k 



2 „ + l_fe_LIi±2| 



2 -!)"((— --)(-!)"- — --) 



2+ir((^-i)(-ir-^-i/4) 



Multiplying a(n) by 4 L — 4 — J , we obtain 1, 3, 20, 272, 7424, . . .. Hence 



1,3,20,272, 



(n+iy 
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k=0 



n 



2k 



1 \ 2n+ l_fc„LH±2j 



fc=0 



77 



2A; 



fc=0 



2 m ^ p +1 ] 



2A- 



That is, the Hankel transform h n (2) of c(n; 2) + c(n + 1; 2) is given by 

2" 



|_2±2j 
k=0 ^ 



2fc 

For our purposes, the following factorization of L _1 is more convenient 
/ 1 ... \ / 1 ... \ / 1 



L" 1 



1 








8 


1 





56^ 


56 


1 


10 


10 


384 


690 


292 


34 


34 


34 







I 

0^0 

i 



/ V 





-8 3 

56 -56 10 

■384 690 -292 34 



/ V 



8 



We now note that the sequence 



fen (2) 

2 m 



is the sequence 62 {n + 1), where b^in) is the 



sequence 1, 3, 10, 34, 116, . . . with generating function 1 _ 4x _|f 2xa and general term 



L-J 

fc=0 



n — k 
k 



c^k ^i—2k ^ y 

k=0 



n — k — 1 

k 



(-2) fc 4 



tin- 2k— 1 



Hence 



h n (2) = 2( n t% 2 (n+ 1). 



Noting that 62 (ft*) is the binomial transform of the Pell A000129(r?, + 1) numbers whose 
generating function is l _^ x _ x2 , we have the following alternative expressions for 6 2 (^) : 



62 (n) 



n / \ k 

e : e 



fc=0 



fc=0 



Eu)E 

i=o 



i=o 
L#J 
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2 2j ~ 



2 fc-2j 



For r = 3, we have 



H 



( 4 15 69 357 

15 69 357 1986 

69 357 1986 11598 

357 1986 11598 70125 

V ; ; ; ; 



\ 



We find that 



L" 1 



/ 1 

\ 

0^0 

^ 

V ; • ; ; 
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198 -131 17 

-2565 2875 -854 73 
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where the sequence 63(71) or 1,4, 17, 73, 314, . . ., A018902 has generating function 1 _ 53 , +33 , 2 
and 



L-J [ n ~ 1 j 

».(») - ±( n ~ k k y-^- 2k -±( n ~ k k ~ i ) 

k=0 ^ ' k=0 ^ ' 



(-3) fc 5 



k r n— 2k— 1 



k=0 



EU)E 

j=0 



k=0 



Eu)E 

j=0 



fc - j 
fc-j 



g2j-fc 



gfc-2j 



Then (n) is the sequence 1, 4, 51, 1971, 228906, .... In other words, we have 

h n {3) = 3( B * 1 )& 3 (n+l). 
We now note that F(2n + 1) has generating function 1 _3~^ 2 with 



-3a;+ 

Ij 

/ n — k\, ^ t „„ oi, fn — k — l 



\knn—2k—l 



F(2n + 1) = £ (-D-3— -SI" : (-1)'3 

fc=0 ^ 7 fc=0 ^ 7 

We can generalize this result as follows. 

Proposition 6. Let h n (r) be the Hankel transform of the sum of the consecutive generalized 
Catalan numbers c(n; r) + c{n + 1; r). Then 

MO = rffl( £ (" " + ') ("OV + 2 ) n " fc+1 - E (" ~ k k ) ("OV + 2)^ 2fc ). 

fc=0 ^ ' k=0 ^ 7 

In other words, h n (r) is the product of r^ 2 ) and the (n + l)-st term of the sequence with 
generating function x^+^+ra 2 ■ Equivalently, 



.n — k+lj ■ / — ' \n — k 

k=0 K 7 k=0 v 

n+l / , , \ k 



\2k—n( \n—k\ 



fc=o v 7 i=o v J/ 

n+l , , ,x L5J 



- 2 



E T Ev 

fc=o v 7 i=o v J 



The two last expressions are a result of the fact that 1 ^( r+2 w+rx 2 * s ^ ne binomial transform 



(r+2)x+rx 2 

1—rx—x 2 



Of , 1 a . 



Berlekamp Massey triangles associated to general- 
ized Catalan numbers 



A natural question that arises when dealing with Hankel matrices is one that is inspired by 
consideration of the Hankel matrix interpretation of the Berlekamp Massey algorithm [2]. 
In our context, this question is that of characterizing the solutions of equations such as the 
following (taking c(n, 3) as an example) 
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\9* J 




I 15 / 



We can define the B-M triangle of a sequence a±, a.2, 0,3, . . . to be the lower triangular matrix 
whose n-th row is the solution of the Berlekamp Massey equations determined by the n-th 
order Hankel matrix of the sequence. 

Example 7. The B-M triangle of the Catalan numbers. We must solve 

(l)(x) = (1) 




and so on. We obtain the triangle 
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-15 7 
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with general term (-l)«-fc((«+ fc + 1 ) - 



) ) and generating function 
1 + x) + xy 



(1 — xy){l + 2x + x 2 — xy) 

Regarding the entries as polynomial coefficients, we see that these polynomials are related 
to the Morgan- Voyce polynomials, themselves a transformation of the Jacobi polynomials. 

In fact, the generating function of the B-M triangle for c(n; r) has generating function 

r(l + x) + xy 



[1 - xy)(l 



l)x + rx 2 — xy) 



These matrices are closely related to the matrices L already studied. The standard Berlekamp 
Massey theory studies the polynomials 



x 



d-i 



i=0 



11 



For the above example, this is 

x 4 - 15x 3 + 75x 2 - U2x + 81. 



Note that the companion matrix of this polynomial is given by 
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In other words, the characteristic polynomial of the last matrix is 



\\ - U2x + 75x 2 - 15x s + x A 



Thus we get 
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3 + .T 



ives 9 — 7x + x l 



ives - 27 + 34x - llx 2 + x 6 



and so on. Forming the matrix of coefficients, we obtain 
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1 



which is the Riordan matrix 



1 



x 



v l + 3x 1 + 4x + 3x 2 ' 

The inverse of this matrix is the L matrix in the LDL T decomposition of the Hankel matrix 
for c(n; 3). The corresponding B-M matrix as defined above is given by the negative of the 
sub-diagonal matrix. 



References 

[1] P. Barry, On integer sequence based constructions of generalized Pascal triangles, 
Preprint, Waterford Institute of Technology, 2005 



12 



[2] N. Ben Atti, G. M. Diaz-Toca and H. Lombardi, The Berlekamp-Massey Algorithm 
Revisited, available electronically at http://hlombardi.free.fr/publis/BMAvar.pdf, 
2006 

[3] A. Cvetkovic, P. Rajkovic and M. Ivkovic, Catalan Numbers, the Hankel Transform, and 
Fiboncci Numbers, Journal of Integer Sequences, 5, May 2002, Art. 02.1.3 

[4] E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, available electronically at 
http : //www . i3s . unice . f r/f psac/FPSAC03/ARTICLES/17 . pdf , 2006 

[5] L. Ferrari, R. Pinzani, Catalan-like numbers and succession numbers, available electron- 
ically at http: //www. arXiv.org/math. CO/0507210,2006. 

[6] J. W. Layman, The Hankel Transform and some of its properties, Journal of Integer 
Sequences, Vol. 4, (2001), Article 01.1.5 

[7] P. Peart, W-J. Woan, Generating Functions via Hankel and Stieltjes Matrices, Journal 
of Integer Sequences, Vol. 3 (2000), Article 00.2.1 

[8] P. Rajkovic, M. D. Petkovic and P. Barry, The Hankel Transform of the 
Sum of Consecutive Generalized Catalan Numbers, available electronically at 
http : //arxiv . org/pdf /math . CO/0604422,2006. 

[9] W-J. Woan, Hankel Matrices and Lattice Paths, Journal of Integer Sequences, Vol. 4 
(2001), Article 01.1.2 



2000 Mathematics Subject Classification: Primary 11B83; Secondary 05E35,11C20,11Y55,15A23. 
Keywords: Hankel transform, Fibonacci numbers, Catalan numbers, Riordan array, Orthog- 
onal polynomials. 



Concerns sequences A000108 , A000129 , A001333 , A007318 , A008288 , A018902 , A079496 



13 



